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Abstract 

Recent experiments suggest that shear-enhanced growth of concentration fluctua- 



tions in polymer solutions is strongly influenced by solid boundaries. We analyze the 
dynamics of a model of a sheared polymer solution, accounting for the effect of stress 
variations on mass flux and for wall-slip. If (and only if) these effects are present, 
the flow exhibits a boundary-localized instability consistent with experimental obser- 
vations. Even in the absence of flow instability, thermally-driven concentration fluctu- 
■ ations are significantly enhanced near the boundary. 

(TJ ! PACS numbers: 61.25.Hq, 83.50.Lh 

The local composition of polymer solutions and blends can be significantly affected by 
flow. In particular, shear-enhanced concentration fluctuations in semidilute polystyrene (PS) 
solutions have been observed via light-scattering and microscopy by a number of groups ([[]; ||; 
El El § El S fD- A P ure ly hydrodynamic mechanism for this phenomenon, in an unbounded 
flow domain, was proposed initially by Helfand and Fredrickson (HF) (|9|) and later extended 
by several authors (10; 11; 12; |l^; [14]). This mechanism is based on the idea that variations 



in polymer stress can drive a flux of polymer molecules - indeed, several authors ( |15| ; |T2 
[T~3| ; |l~6| ) have developed simple, two-fluid models which demonstrate quite generally that the 
polymer flux in a solution, j, has the general form j = -DVn + (Z)/&bT)V • r, where n is 
the polymer number density, r is the polymer extra stress tensor, D is the polymer gradient 
diffusivity, ks is the Boltzmann constant, and T is the temperature. Furthermore, a detailed 
kinetic theory derivation for a dilute solution ( |P7|) leads to a very similar expression for the 
flux. A scattering peak is predicted to arise when the diffusion and relaxation times are 
equal, which occurs at a wavelength of 2ir/k m \J DX which is, not coincidentally, the only 
inherent length scale in the problem. 

This length scale is observed in light scattering experiments. Hashimoto and coworkers ([!]; 
H]; II H) have performed a series of experiments with PS dissolved in dioctylphthalate (DOP) 
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and have measured the scattering intensity by illuminating along the shear-gradient direction 
(the ^/-direction of Fig. [I]). They observed a "butterfly" scattering pattern indicative of 
flow structures with wave vectors roughly aligned in the flow direction. The wavenumbers 
obtained from the patterns were on the order of k ~ 0.6/im -1 giving a length scale for the 
fluctuations of approximately 10/im. Estimating the relaxation time as the reciprocal of the 
shear rate for the onset of shear thinning (§) and using the diffusivity measured by Wu et 
al. (HI) for a similar system gives a predicted wavenumber of k « 0.7/mi _1 , which is close to 
the measured value. This length scale was also observed by Wu et al. (Q), who measured the 
scattering intensity by illuminating along the vorticity direction (the z-direction in Fig. H). At 
low shear rates, the scattering peak was located at k ~ 10/im -1 , which is in good agreement 
with their predicted value of 9.4//m _1 . 

Even though both of these experiments show peaks at the requisite length scale, quali- 
tative differences exist in the shear rate dependence of the results. In particular, Hashimoto 
and coworkers, who measured the intensity essentially averaged over a volume containing 
both near-surface regions, found that the intensity increased markedly above a critical shear 
rate, which was roughly the reciprocal of the polymer relaxation time (an approximately 100- 
fold increase in scattering for a 10-fold increase in shear rate). Conversely, Wu et al, whose 
measurements sample only the bulk, not the near-surface regions, found a much smaller 
increase in intensity as the shear rate was increased, approximately an 11-fold increase in 
scattering intensity for a 25-fold increase in shear rate. 

A recent experimental study by Mhetar and Archer ( Jl8| ) has demonstrated several novel 
features of shear enhanced concentration fluctuations in entangled polymer solutions. They 
studied semidilute solutions of PS in diethylphthalate (DEP), to which tracer particles were 
added to measure the velocity profile, using a planar Couette cell with (high energy) silica 
walls. In this situation, wall-slip occurs primarily due to the disentanglement of bulk chains 
from those adsorbed to the walls, rather than to the formation of a depletion layer, as arises 



in dilute solutions (IS; EOT). Two features of the flow of these solutions stand out. First, no 



turbidity was observed at low shear stresses where slip was weak; enhanced concentration 
fluctuations were observed, however, at large shear stresses, where slip was significant. The 
length scale for the fluctuations was observed to be ~ 10/im, as was the extrapolation, or 
slip, length measured for these solutions. Second, the concentration fluctuation enhancement 
began near the boundaries and modification of the surface to increase slip delayed the onset 
on fluctuations to much larger shear stresses. These observations strongly suggest that slip 
plays a role in the formation and development of these fluctuations. 

Taken together, experiments suggest that the HF mechanism, while being able to explain 
bulk behavior, only partially accounts for the physics leading to enhancement of fluctuations 
near surfaces, a conclusion that can be put into perspective by consideration of the fluid 
dynamics of polymer flows near boundaries. Gorodtsov and Leonov (^) analyzed the plane 
Couette flow of an upper convected Maxwell (UCM) fluid (a model of a melt) and found that, 
while the flow is always linearly stable, the most slowly decaying modes are localized near 
the bounding surfaces. As a result, we might expect the formation of boundary localized 
structures due to random concentration fluctuations in the polymer solution. In addition, slip 
has been shown to lead to boundary-localized flow instability in polymeric simple shear flows 
under certain conditions (|23| ; |24] ; [25]) . This suggests that boundary-localized concentration 
fluctuation dynamics can appear by analogous routes. Our goal here is to analyze the 
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coupling between concentration, stress, and slip in bounded shear flow, to elucidate the 
effects of boundaries on the dynamics. 

We consider flow in the plane Couette geometry shown in Fig. [l|. To describe the behavior 
of the solution, we use the kinetic theory model derived by Beris and Mavrantzas (|26|), for a 
solution of Hookean dumbbells, and written in terms of the conformation tensor, a = (QQ), 
which is the second moment of the orientation vector of the dumbbells, Q. Nondimension- 
alizing stresses with n^ksT, where Uq is the average concentration of the solution, lengths 
with the characteristic length \J DX, velocities with 7*\/.DA, time with A, and concentration 
with rig gives the following model equations: 

Oct 

— + v ■ Vck - cx ■ Vv - (Vv) T -ol + a-nd = V 2 a, (1) 

— + v ■ Vn = -V ■ (-Vn + V ■ r) , (2) 

-Vp+V -t + SV 2 v = 0, (3) 
V • v = 0, (4) 

where v is the velocity, S is the ratio of the solvent to polymer viscosities, n is the number 
density of dumbbells, r = a — nS is the polymer extra stress and p is the pressure. In 
presenting the model, we have suppressed any dependencies of properties on concentration - 
in the analysis below we only consider small perturbations from spatially uniform stress and 
concentration fields, in which case the terms that reflect these dependencies do not appear. 
Eq. ffl describes the evolution of polymer conformation and is a generalization of the upper 



convected Maxwell (or Oldroyd-B) model (27) to include spatial diffusion of conformation 



(on the RHS). In most cases, this term is negligible (|2"6|), but here, the flow structures we are 
interested in have a length scale set by the diffusivity. Eq. ^| is the conservation equation for 
the dumbbells, explicitly showing the dependence of the polymer flux on stress variations. If 
diffusion of the polymer molecules is neglected, the right-hand sides of Eqs. [I] and [| vanish, 
recovering the Maxwell model. Eqs. |^ and [| describe conservation of momentum and mass. 
Other phenomenologically motivated models for such systems are virtually identical to this 
one, but do not generally include the conformational diffusion term (|]; |10[ |12; [13|; [14]) , 



the primary difference being the use of the polymer phase velocity rather than the mass- 
average velocity in the deformation terms of the equation for polymer conformation - our 
preliminary results with such a full two-fluid model yield qualitatively identical results to 
those presented below. Furthermore, all earlier works have restricted themselves to bulk 
behavior and We < 1. 

Because of the high energy surfaces in Mhetar and Archer's experiments, it is reasonable 
to assume that the adsorbed chains do not desorb to a significant degree in flow, and that 
no additional bulk chains adsorb. Therefore, we apply a no-flux boundary condition for the 
polymer concentration: n ■ j = dn/dy — dr yx /dx — dT yy /dy — dr zy /dz = 0, where n is the 
unit normal on the boundary. Inclusion of the diffusion term in the constitutive equation 
necessitates the specification of boundary conditions for the conformation tensor. For the 
same reason as above, we assume that there are no sources or sinks for polymer conformation 
at the walls, giving a no-flux boundary condition for the conformation, da/dy = 0. We have 
also performed computations using the classical UCM equation as the boundary condition 
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for ol (|2q) and found no qualitative difference in the results. 

The model is completed by specifying the slip boundary condition at the interface, which 
crudely reflects the interaction between the bulk chains and those attached to the wall. We 
use here the simplest nontrivial model, the Navier slip boundary condition u s = er yx (= b = 
b* /V DX at steady state); e is given by bjWe. 

We have studied the dynamics of this model primarily by linear stability analysis of the 
plane Couette base state, given by a = {v x ,v y , v z , f xx , f yx , f zx , f yy , f zy , f zz , n, p) 
= (W(y + b), 0, 0, 2W 2 , We, 0, 0, 0, 0, 1, 1), where ex. has been eliminated in favor of r. Small 
perturbations are added to this base solution, i. e. a = a(y) + a(x,y, z,t), and only the 
terms linear in the perturbation are retained in the model equations. The perturbations 
take the normal mode form: a(x,y,z,t) = a(y)e th ' x e~ lat + c.c, where k = (k x ,0,k z ). This 
formulation yields a generalized eigenvalue problem for the eigenvalues {a} with associated 
eigenfunctions {a(y)}. If Im(o~) > then infinitesimal disturbances grow and the flow is 
unstable. 

Prior work (24; 25) with the UCM model and preliminary work with the present model 
have shown that the eigenfunctions of the discrete spectrum are localized near the boundaries. 
Furthermore, in experiments, \J DX <C I* . Therefore, we solve the eigenvalue problem on a 
semi-infinite domain (SID). We map y G [0, oo) onto the computational domain, £ G [—1, 1], 
using the mapping ?/ = 0.1(l + £)/(l — £), and then perform Chebyshev collocation on the 
mapped system (|29"D , using N + 1 collocation points. The eigenvalue problem is solved using 
a public domain routine (|30|). 

To begin the discussion of the results, we note that if the possibility of concentration 
fluctuations is suppressed in the model (by setting D to zero, thus eliminating the RHS 
of eqs. |l] and |2|), plane shear flow with either the no-slip ( p2]) or Navier slip ( PH |24|; ^5|) 



boundary condition is stable for all Weissenberg numbers. In addition, we have found that 
if concentration variations are permitted but slip is not (by setting 6 = 0), the flow is also 
stable. In contrast, if slip and concentration variations are allowed, the flow can become 
unstable when We exceeds a critical value We c - a new class of polymer flow instabilities 
arises. Fig. || shows neutral curves, plotted as We c vs. k x for two dimensional disturbances, 
i.e. k z = 0, with b fixed. For b = 0(1), the critical wavenumber is 0(1), so that the critical 
length scale for the instability is V DX. We have also examined stability with respect to three 
dimensional disturbances ik z ^ 0), but find that two dimensional disturbances are always 
the most dangerous. At low wavenumbers (those most easily observed by microscopy), We c 
first decreases and then increases with increasing b - the increasing section of the curve is 
in agreement with the observations of Mhetar and Archer (|1^) that treating the surface to 
increase slip delayed the onset of fluctuations. Note that there is a transition from one mode 
of instability to another, leading to multiple minima in the neutral curve, as most clearly 
seen in the curve for b = 2. For larger values of b, the neutral curves collapse toward a single 
curve. Similar results are found for larger values of S*, but the critical Weissenberg numbers 
are much lower, We c ~ 2 for S > 1. (However, for S = 0(1), the solution would typically be 
too dilute for slip to be a realistic boundary condition.) 

Fig. |3] shows a density plot of a typical destabilizing concentration perturbation. Near 
the boundary, the local wavevector has the same orientation as that predicted by the HF 
mechanism (|| [13]). In agreement with the observations of Mhetar and Archer (|T8|), pertur- 
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bations are localized near the surfaces and are observed at large We. As a result, we would 
expect scattering to be more pronounced in the boundary regions than in the bulk. 

We have also investigated the dynamics of Brownian fluctuations in this model. Ji and 
Helfand ( |l3l) have demonstrated that the velocity is a fast variable and therefore, since 
introduction of concentration fluctuations into this system results in adding fluctuations 
directly to the stress, we only incorporate the random forcing into the polymer mass balance. 
This is carried out by adding a term V ■ w to Eq. ||], where w is the random contribution 
to the polymer flux and satisfies (w(x, t)w(x', t')) = [(2n*(DX) 3 / 2 )]6 S(x - x')5(t - if) 
To conserve mass, we require w y = on y = 0, 1. We follow the common procedure (|| 



10| ; [12| ; p~3| ) of assuming the noise is small and contributes only to the linearized equations. 
Again, Chebyshev collocation is used for the spatial discretization in the y-direction with a 
Fourier transform in the x-direction, but because the collocation points are not uniformly 
distributed, the random forcing is weighted so that it is uniform in magnitude throughout 
the domain. The resulting system of equations was time-integrated on a bounded domain 
using an implicit Euler scheme with the top plate located at I* = 100V-DA. Results are 
collected only after the initial transient has decayed and fluctuation effects dominate. 

To quantify the response as a function of time, we compute the concentration correlation 
function for a given x-wavenumber k x : (n(x; k x )n(x'; k x )) = 2Re(h(y; k x )n*(y'; k x )) cos{k x {x— 
x')) — 2lm{h(y; k x )n*(y'; k x )) sm{k x (x — x')). Fluctuations are strongly correlated near the 
surfaces, even at equilibrium {We = 0). At equilibrium, lm(n(y)n* (y')) is zero, but for 
nonzero We, it assumes values consistent with enhanced fluctuations whose wave vectors are 
rotated from the flow direction into the first quadrant, while the magnitude of the real part 
remains essentially unchanged. A snapshot of a typical concentration profile is shown in 
Fig. |] for W = l,k x = l,b = 1 and the orientation, which is again consistent with the bulk 
HF mechanism, is clear. Hence, fluctuations form near the surface by mechanisms analogous 
to those in the bulk, but are enhanced because of the vicinity to the no-flux boundary. Re- 
sults for the no-slip case are virtually identical - the response of the flow to Brownian noise 
is quite insensitive to the presence of slip. Therefore, even in the absence of flow instability, 
(i.e. when We < 0(1) and/or slip is absent), the near surface regions may make a nontrivial 
contribution to the scattering signal, particularly its anisotropy under flow. The insensitivity 
to slip is consistent with what is known about the eigenspectra of viscoelastic flows: whether 
or not there is slip, the slowest decaying eigenmodes are localized near solid boundaries. 

These results clearly suggest that the bounding surfaces play an important role in the 
formation of enhanced concentration fluctuations in sheared polymer solutions. In the con- 
text of the present model, slip leads to a new class of viscoelastic flow instabilities which 
result from the interaction of slip with stress and concentration, and even in the absence 
of slip, Brownian fluctuations are selectively enhanced near the surfaces. These results are 
consistent with experimental observations and underscore two points regarding the flow be- 
havior of polymeric liquids: (1) the distinctness and importance of the dynamics of flowing 
polymers near boundaries, even at the continuum level, and (2) the importance of couplings 
between various phenomena for the dynamics and stability of these flows. 

The authors gratefully acknowledge the National Science Foundation and the Nontenured 
faculty grant program of the 3M Company for supporting this work. 
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Figure 1: Plane Couette geometry showing slip between the solution and solid surface. The 
flow is in the x-direction, y is the gradient direction and z is the vorticity or neutral direction. 
I* is the gap width, u* is the slip velocity, 7* is the true shear rate, and b* = u*/Y is the 
steady-state extrapolation length. The asterisks (here and in the text) denote dimensional 
quantities. 




Figure 2: Neutral curves for fc 2 = 0, 5 = 10 2 ,iV = 96, and various values of the extrapola- 
tion length. O - b = 1, Q - b = 2, □ - b = 4, A - b = 8, and v - b = 10. 
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Figure 3: Concentration component of the unstable eigenfunction; k x = 0.4, b — 10, We — 10, 
S = 0.01, JV = 96. 
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Figure 4: Typical snapshot of a concentration profile driven by random fluctuations. The 
parameters are: We = 1, k x = 1, b = 1, S = 0.01, iV = 192. 
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